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O I Measured values of the brightness temperature of low-frequency synchrotron radiation emit- 
^ ■ ted by powerful extragalactic sources reach 10^^-10^^ K. If some amount of nonrelativistic 
;_i . ionized gas is present within such sources, it should be heated as a result of induced Compton 
^ scattering of the radiation. If this heating is counteracted by cooling due to inverse Compton 

scattering of the same radio radiation, then the plasma can be heated up to mildly relativis- 
I tic temperatures kT ~ 10-100 keV. The stationary electron velocity distribution can be 
either relativistic Maxwellian or quasi-Maxwellian (with the high- velocity tail suppressed), 
depending on the efficiency of Coulomb collisions and other relaxation processes. We derive 
several easy-to-use approximate expressions for the induced Compton heating rate of mildly 
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Tj- ! relativistic electrons in an isotropic radiation field, as well as for the stationary distribution 

^ I function and temperature of electrons. 

^ ' We give analytic expressions for the kernel of the integral kinetic equation (one as a function 
of the scattering angle and another for the case of an isotropic radiation field) , which describes 

^ . the redistribution of photons in frequency caused by induced Compton scattering in thermal 

^! plasma. These expressions can be used in the parameter range hu <^ kT <^ O.lmc^ (the 

O I formulae earlier published in Sazonov, Sunyaev, 2000 are less accurate). 



1 Introduction 

Milliarcsecond-resolution interferometric radio observations of the central parsec-scale re- 
gions of active galactic nuclei often reveal sub-structure emitting low-frequency radiation 
with inferred rest-frame brightness temperatures Th ~ 10^^-10^^ K. If there were thermal 
plasma present in these regions, it could be heated efficiently as a result of induced Compton 
scattering of the radio radiation on the electrons (Levich, Sunyaev, 1971). 

Such a plasma has not been observed so far in extragalactic sources. On the other hand, 
there is evidence that relatively cold matter constitutes a significant fraction of the total mass 
contained in the jets, including their inner (< 10^^ cm) regions, of the famous galactic source 
SS433 (see Vermeulen, 1992 for a review). Of course, these jets are only mildly relativistic 
[v ~ 0.26c) and therefore quite different from the relativistic jets in active galactic nuclei 
(which have bulk Lorentz factors 7 ~ 5), but it is reasonable to expect that some amount 
of quasi-thermal plasma may be present in extragalactic jets as well (see, e.g., Celotti et al., 
1998). 



One can derive directly from the Kompaneets kinetic equation the induced Compton heating 
rate of thermal electrons located in an isotropic radiation field (Levich, Sunyaev, 1971). The 



resulting expression is, however, only applicable when the electron temperature is relatively 
low, kT ^ a few keV. As the electrons become more relativistic, their heating rate monoton- 
ically decreases. Theoretical efforts in the 70-s (Vinogradov, Pustovalov, 1972; Blandford, 
1973; Blandford, Scharlemann, 1975) resulted in a number of useful formulae that allow 
one to determine the heating rate of relativistic electrons and their distribution function in 
various limits, e.g., for ultrarelativistic electrons, for narrow radiation beams, or for some 
specific radiation spectra. Illarionov and Kompaneets (1976) have derived a general expres- 
sion that gives the heating rate for an electron moving with an arbitrary velocity in a given 
isotropic radiation field. This formula is, however, rather complex: one needs to compute a 
3-dimensional integral in order to find the heating rate of an ensemble of electrons with a 
given velocity distribution (e.g, Maxwellian) . 

In the present paper we demonstrate that electrons can be heated by Compton scattering 
up to mildly relativistic temperatures kT ~ a few tens of keV, but not more, in an isotropic 
radiation field with Tj, < 10^^ K. For these plasma temperatures, the nonrelativistic esti- 
mates for the heating rate and some relevant quantities become inaccurate. However, these 
expressions can be modified by adding to them a few correction terms, thereby retaining the 
original simple structure, as shown below. 

It should be noted that although we mentioned above only active galactic nuclei, the induced 
Compton effect may play a major role also in other astrophysical situations. We also note 
that the problem of plasma heating near, but outside, a source of low-frequency radiation 
requires a special study. 

2 Heating and cooling of thermal electrons during Comp- 
ton scattering 

Let us derive the rates of heating and cooling of mildly relativistic thermal electrons {kT ^ 
O.lmc^, where T is the electron temperature) in an isotropic radiation field, as a result of 
both the spontaneous and induced Compton scattering processes. 

The radiation field may be defined by its spectral energy density, e^, (measured in units 
of erg cm~^ Hz~^), or, equivalently, by the occupation number in photon phase space, 
n — ei,c^ /{STihu^). We require that the spectrum be sufficiently broad, which means 
that its characteristic width Au/i/ should be much larger than the typical fractional fre- 
quency shift acquired by a photon during a scattering event either due to the Doppler effect, 
~ ±(kT/mc^y^'^, or due to Compton recoil, ~ —hv/mc^. In this case the integral kinetic 
equation describing the time evolution of ej, (caused by the Compton interaction of the radi- 
ation with the electrons) can be transformed into a Fokker-Planck-type differential equation. 
In the limit of kT, hv <^ mc^, the resulting equation is the Kompaneets (1957) equation. 

From the Kompaneets equation one can derive the expressions for the electron heating and 
cooling rates (Levich, Sunyaev, 1971), which are fairly accurate for electron and photon en- 
ergies ^ O.Olmc^. The Kompaneets equation can be extended into the mildly-relativistic do- 
main {kT, hv <^ O.lmc^) by adding to it relativistic correction terms ~ {kT /mc^)'^{h'u/m(?)'^ 
(Itoh et al., 1998; Challinor, Lasenby, 1998; Sazonov, Sunyaev, 2000). For example, the 
first-order generalization of the Kompaneets equation (which obtains by Fokker-Planck-type 
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expansion of the integral kinetic equation to the fourth order in Au) is 
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where ctt is the Thomson cross section and A^e is the number density of electrons. 

Multiplying Eq. (|I|) by hv and integrating over the photon phase space leads to an expression 
for the variation of the radiation total energy density, e = /q°° e^dv, and, at the same time, 
for the change of the mean electron energy, E^, ^ 1.5kT{l + 1.25kT/mc'^): 
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The first two terms in the square brackets in Eq. (§) describe the electron heating and cooling 
rates in the nonrelativistic limit (Peyraud, 1968; Zel'dovich, Levich, 1970; Levich, Sunyaev, 
1971). The additional terms represent relativistic corrections of the order of kT/mc^ or 
hu/mc^. Note that the correction to the heating rate associated with the induced Compton 
effect (the corresponding terms are nonlinear in n) includes a term which is expressed through 
the square of the derivative, [dnjdvY- Higher-order [{kT/mc^Y, etc.] corrections to the 
induced- Compton energy transfer rate, which are not given in Eq. (Q), depend on higher- 
order derivatives of n, as is explicitly shown below (in §3). 

Using Eq. (^, we can, for example, find the energy transfer rate between electrons and 
black-body radiation of temperature T^, 

de , , . / kT kTr \ ( 5 kT , fcT^ 
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dt \mc^ mc' J \ 2mc'' mc' 

This expression was originally derived (in a different way) by Woodward (1970). 

We note that the problem of the energy transfer between radiation and a thermal plasma 
via spontaneous Compton scattering is well studied, and the corresponding part of Eq. (|) 
follows directly from the first moment of the Compton scattering kernel (Shestakov et al., 
1988; Sazonov, Sunyaev, 2000; also see the Appendix). We are concerned here with the 
terms describing the contribution of induced scattering. 
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3 Induced Compton heating of thermal electrons 



Powerful extragalactic radio sources emit low-frequency continuum radiation that in some 
cases has a very high brightness temperature, Tb = nhu/k ~ 10^^-10^^ K, so that kTb > mc^. 
Induced Compton scattering can be a major mechanism of heating of free electrons situated 
in such a radiation field. Moreover, nonrelativistic estimates lead to the conclusion that 
electrons can be heated up to relativistic temperatures kT ~ a fraction of kTf, (Levich, 
Sunyaev, 1971). Clearly, nonrelativistic treatment is inappropriate here. Below we show 
that in reality, electron temperatures achievable in this situation are only mildly relativistic 
(i.e., kT is significantly less than mc^). As a result, all relevant physical quantities can be 
described by simple analytic expressions. 

We consider throughout an isotropic radiation field. In this paragraph we also assume a 
relativistic Maxwellian distribution of electrons (which can be maintained, e.g., by Coulomb 
collisions). This last assumption is dropped in §4. 

The rate at which energy is transfered by induced Compton scattering from an isotropic radi- 
ation field to a single electron moving at a speed v = (3c\s given by (Illarionov, Kompaneets, 
1976) 
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Eq. (^) is valid for arbitrary electron velocities, provided that hu <^ mc^. This formula can 
be simplified in some limits, in particular, for spectrally narrow radiation lines (Az/ <C v), 
such as those produced in cosmic masers and laboratory lasers (Vinogradov, Pustovalov, 
1972; Blandford, Scharlemann, 1975). We are instead interested in the Au ^ v case. 

In order to calculate the instantaneous heating rate of an actual plasma, it is necessary 
to weight the result of Eq. (^) with a given distribution of electron velocities, f^. The 
computation procedure to determine the heating rate will therefore be three-dimentional 
integration: 1) over the radiation spectrum [when calculating $(?/')]! 2) over d^' and 3) over 
However, if the electrons obey a relativistic Maxwellian distribution and are only mildly 
relativistic [kT ^ O.lmc^), it is possible to considerably simplify the computation procedure, 
by reducing it to integration over the radiation spectrum, as demonstrated below. 

When kT <^ O.lmc^, the majority of electrons have (3 <^ 0.5, and therefore the typical random 
Doppler frequency shift caused by a scattering is relatively small: ^ u. The following 



4 



Taylor expansion of the function entering Eq. (^) is then justified: 
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A similar expansion is possible for the kernel G{(3, (3') of Eq. 
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The subsequent straightforward calculation leads to 
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Note that the Planck constant, h, is absent from Eq. which reflects the well-known fact 
that the process of induced Compton scattering is classical (see, e.g., Zel'dovich, 1975). 

Our next step is to replace in Eq. (H) by their average values {P'^), {P^), 

calculated for a mildly-relativistic Maxwellian distribution at temperature T, which is given 
by (e.g., Sazonov, Sunyaev, 2000) 
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The resulting heating rate as a function of the plasma temperature is 
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where Jo, Ji, J2 and J3 were introduced in Eq. (H). 



In the case of cold electrons [kT <^ mc^), only the leading term of the series in powers of 
kT/mc^ in Eq. (|10D is important, and the heating rate due to induced Compton scattering 
is described by the well-known (Zel'dovich, Levich, 1970; Levich, Sunyaev, 1971) formula 
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Both the nonrelativistic expression ( [TlD and the first-order relativistic correction to it (see 
Eq. []10[) were already obtained in §2 following a different approach; see Eq. (^. Thus, two 



independent methods give the same result. 



3.1 Heating in a field of self-absorbed synchrotron radiation 

Let us now consider one particular spectral distribution, namely the spectrum of synchrotron 
radio emission with self-absorption at low frequencies. The radiation spectrum generated by 
a spherically-symmetric homogeneous source is (Gould, 1979) 
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The shape of the distribution (1121) is defined by a single parameter, p, which is the index 
of the power-law energy distribution of the electrons producing the synchrotron radiation, 
dNe ~ j'^dj. The other two parameters that appear in Eq. ([I^), and A, determine 
the position of the spectrum along the frequency axis and its amplitude, respectively. Far 
enough from the peak of intensity (fpeak ~ l-15z/o), the spectrum (0) assumes a power-law 
shape: e[J ~ z/^/^ in the region of self-absorption <^ Uq) and e° ~ iiC^-p)/"^ in the optically 
thin part {u ^ Uq). 

Since the integral e = e^du diverges at cxo for p < 3, and it is this quantity that determines 
the Compton cooling rate, which will be taken into account below, we modify Eq. (|T2p as 
follows: 

e.(-) = | ?^r^\ 5 0, , (13) 

Here we have assumed that a steepening of the spectrum (increase in the slope by 0.5) takes 
place above some "break" frequency, Ub ^ z^o, due to the fast synchrotron cooling of more 
energetic electrons. 

To be concrete, we take the slope of the optically thin part of the spectrum (prior to the slope 
break) to be a = —0.7, which corresponds to p = 2.4. We additionally adopt Ub = lO^t'o as a 
fiducial value for our treatment of electron cooling, although this parameter just enters as a 
multiplicative factor in the relevant formulae as long as Ub ^ i^o- The resulting spectrum is 
plotted in Fig. |I|. Note that the low-frequency breaks observed in the spectra of some radio 
sources are sometimes interpreted as being caused by mechanisms other then synchrotron 
self-absorption, which include free-free absorption in the ambient medium (e.g., Bicknell et 
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al., 1997) and induced Compton scattering either outside or inside the radio source (e.g., 
Sunyaev, 1971; Sincell, Krohk, 1994; Kuncic et al., 1998). However, the radiation field inside 
the source, which we are concerned with, may well be the superposition of self-absorbed 
synchrotron spectra generated by individual plasma components. 

Fig. 1^ shows the heating rate of thermal electrons exposed to self-absorbed synchrotron 
radiation as a function of the electron temperature. The exact result was obtained by 
weighting the heating rate for monoenergetic electrons, given by Eq. (P, with the rela- 
tivistic Maxwellian distribution, dN^ = const 7(7^ — 1)^^^ exp (— 7mc^//i;T)(i7. We have 
additionally verified this dependence with Monte-Carlo simulations (using the Comptoniza- 
tion code described in Pozdnyakov et al., 1983 with a slight modification to allow for induced 
Compton scattering). The computation proves to be faster with the semi-analytic formula 
(P of lUarionov and Kompaneets. Also presented in Fig. ^ are various approximations 
for the heating rate that result from retaining a different number of temperature terms in 
Eq. (0). 

One can see that the heating rate begins to decrease appreciably (> 5%) at A;T ~ 5 keV. The 
temperature domain where this decrement is described very well by the approximate relation 
( p!OD extends to kT ~ 30 keV, when the heating rate is already smaller by ~ 30% than in the 
case of cold electrons. The further decrease in the heating rate that takes place at yet higher 
temperatures cannot be described properly in the framework of the Fokker-Planck approach, 
which led to Eq. ([T0|) . We found it convenient to describe the exact dependence presented 
in Fig. ^ by an approximate formula, which is accurate to within 3% for kT < bmc^: 
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where, as follows from Eq. ([1^), ( p!^ : 



= 1.1 ■ lO^MVo-^eVs"^ 



(15) 



{A is measured in erg Hz ^ cm ^ and vq in GHz). 

It is worth mentioning that the heating rate approaches an asymptote W^{T) ~ {mc^/kT)^ 
when kT ^ mc^. This results from the fact that only the lower-energy part of the relativistic 
Maxwellian distribution, i.e., electrons with 7^1, significantly contributes to the net heating 
rate (because ~ 7"^ for 7^1 — see Illarionov, Kompaneets, 1976), and the relative 
number of such electrons is proportional to (mc^/A;T)^. 

Although we have assumed a particular slope (—0.7) for the optically-thin part of the radi- 
ation spectrum, it turns out that the dependence of the heating rate on the electron tem- 
perature changes very slowly as the slope varies. Quantitatively, W~^{T) remains the same 
to within 10% for a in the range [—0.9; —0.5]. Thus, formula ( [T^ ) is quite useful in that it 
allows obtaining reasonably good estimates for the heating rate of substantially relativistic 
Maxwellian electrons in an isotropic field of self-absorbed synchrotron radiation. 

Of course, the mildly-relativistic Eq. (|TUp, which is applicable when kT <^ 30 keV, explicitly 
depends on the spectral distribution, and so can be used to calculate the heating rate for an 
arbitrary (broad) radiation spectrum. 
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3.2 Stationary temperature of electrons 



Zel'dovich and Levich (1970) studied, in the nonrelativistic approximation, the problem 
about the estabhshment of a stationary distribution of electrons during their interaction 
with an isotropic field of low-frequency radiation of high brightness temperature. These 
authors showed that if the induced Compton heating is counteracted by cooling due to inverse 
(spontaneous) Compton scattering, then the stationary distribution will be Maxwellian with 
a temperature 

where e = /q°° ey dv. 

It is reasonable to suggest that in the mildly-relativistic regime there will be only small 
deviations from a Maxwellian distribution. Furthermore, if other mechanisms, e.g.. Coulomb 
collisions should play a significant role in the redistribution of electrons in momentum space, 
then a Maxwellian distribution can be maintained easily. We postpone a detailed discussion 
of questions related to the electron distribution until §4. Here we will continue to assume 
as before that the electrons obey a relativistic Maxwellian distribution, and will find the 
stationary electron temperature, considering inverse Compton scattering (of the same low- 
frequency synchrotron radiation that is heating the electrons) the only cooling agent. Other 
possible cooling mechanisms, which may prove more important under certain conditions, will 
be mentioned in §3.3. 

The inverse Compton cooling rate is given by (e.g., Pozdnyakov et al., 1983) 

W^-m=((7) + ;^)w^o-(r), (17) 

where the mean electron energy (in units of mc^) is 

^kTK2{mc^/kT) + mc^Ki{mc^ /kT) ^, 



(7) 



2kTKi{mc^/kT) + mc^ Ko{mc^ / kT) ^ 



Kp{x) are modified Bessel functions, and the cooling rate in the nonrelativistic limit {kT ^ 
mc^) is 

W,;(T) = 1^. (19) 
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We can expand Eq. (0) in powers of kT/mc^ to obtain a formula applicable in the mildly- 
relativistic limit, which is similar in structure to the corresponding relation for the heating 
rate (Eq. [H 
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An excellent fit to the exact formula ([17|) in the range kT ^ 5mc is provided by 
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with Wq given by Eq. (pj|). 

The dependences W~{T) described by Eqs. (p!7D and (^) are plotted in Fig. ^ and can be 
compared to the corresponding dependences for the heating rate (Fig. ||a). We note that 
the relationship W~{T), of course, does not depend on the shape of the radiation spectral 
distribution (depending only on the total energy density), as opposed to W~^{T). 

We can now find the equilibrium electron temperature by solving the equation W~^{Teq) = 
W~{Teq). In order to proceed, we need to specify the constant A appearing in Eq. (|T2]), which 
defines the amplitude of the spectral distribution. It is natural to express this coefficient 
through the maximal brightness temperature, T™^^ = max[e^{i>)c^/{87ri''^)]. peaks at 
j^max _ g g]^jyp fQj^ Q, — —0.7 (although the position of the peak is only marginally dependent 
on a for typically observed radiation spectra of extragalactic radiosources) . We find: 

A ^ 22-nkT^''''vl/c^. (22) 



Fig. ^ shows the equilibrium electron temperature, Teq, as a function of T™^^. The exact 
result is compared with the nonrelativistic result (Eq. |jT6|) and different-order mildly- 
relativistic estimates that were obtained by equating Eq. (^) and Eq. (pOl). One can see 
that the nonrelativistic Eq. (0) is still valid for Teq ^ 5 keV, which corresponds in our case 
to brightess temperatures T™'^^ ^ 3 ■ 10^° K. In this regime 

T-imax / , . \ —0.33 

fcT° = 1.9-^ ( ) keV. (23) 
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The next decade of values of the equilibrium temperature, up to kT^^ ~ 40 keV, is well 
described by the approximate formulae (p!0D and (^0]). Note that fcTeq = 40 keV corresponds 



to Tfo"^^^ ~ 4 ■ 10"(i/b/1000i/o)°-^^ K. An important conclusion can be made: electrons can be 
heated up to mildly relativistic temperatures, kT ~ a few tens of keV, as a result of induced 
Compton scattering of synchrotron radiation with Tb ~ 10^^-10^^ K, but not above these 
temperatures. 



3.3 Evolution of the electron temperature during Compton inter- 
action 

Let us now address a related timing problem, namely, examine how rapidly electrons can be 
heated to mildly relativistic temperatures through induced Compton scattering. 

We will first estimate the basic characteristic quantities in the nonrelativistic limit. Substi- 
tuting A given by Eq. ( P^ into Eq. (|r^), we find 

HV^1.4.10-»(^) (^)eVs-'. (24) 

If there were no cooling, then initially cold electrons would acquire a kinetic energy of kT^^, 
given by Eq. (p3D , during 

/ T^max N -1 / \ -3 / \ -0.33 

s. (25) 
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We notice that the heating time depends strongly on the characteristic frequency of the 
synchrotron self- absorption: theat ~ t'o^^- 

Using Eq. (|19D, we can find the corresponding cooling rate: 

In the mildly-relativistic regime, the heating and cooling rates will be respectively smaller 
and larger than those given by Eqs. ( p4|) and (^61) . We have computed the evolution of 
the electron temperature for a set of T™^^ values by integrating the following differential 
equation: 

dT _/ d{j){T) y' W-^{T)-W-{T) ^ 
dt \ dT J mc^ ' 

using Eq. ( pISD to represent the dependence (7) (T) and the fitting formulae ( p!4D and (|2T|) for 
W~^{T) and W~{T), respectively. The resulting time histories are presented in Fig. ^, for 
which the value z/q = 1 GHz was used. For a given value of z/q, one should simply rescale the 
time axis in Fig. ^ as (z/q/IGHz)"'^. 

Each of the time histories presented in Fig. ^ can be divided into two intervals. During the 
earlier period, the temperature grows linearly because ^ W~ . As T becomes ^ O.STeq, 
the second (longer) period starts, during which the cooling plays a major role and the 
temperature slowly approaches the equilibrium value. This transition period is additionally 
lengthened by relativistic effects (compare the different solutions in Fig. ^). We can define 
a characteristic time of induced heating as the time needed to heat the plasma to ~ O.STeq. 
To give an example, for T^^^ = 10^^ K, a time of ~ 3 ■ 10^(z/o/lGHz)~^ years is required 
to heat the electrons to kT = O.S/cTcq ~ 34 keV. If further, uq ~ lOGHz, the heating time 
becomes ~ 30 years. Interestingly enough, the simple nonrelativistic Eq. (^5]) provides a 
good estimate (within a factor of 2) for the heating time even for values of T™^^ as high 
as ~ 10^'^(z/f,/1000z/o)°'^'^ K (of course, the equilibrium temperature in this case is much less 
than the nonrelativistic estimate). 

In real situations, there may be mechanisms operating by which the plasma cools more ef- 
ficiently than by inverse Compton scattering of the synchrotron radiation. One should then 
modify accordingly the cooling rate W~{T) in our treatment above. If the energy density 
of a possible high-frequency radiation component is larger than that of the low-frequency 
synchrotron emission, then the contribution of this component to the inverse Compton cool- 
ing rate will be accordingly larger. Also, cooling due to free-free transitions will become 
more important than inverse Compton cooling if the plasma is dense enough, namely when 
Nf.T~^^'^e~^ > 10^ K~^/^ erg~^. For the synchrotron spectrum described by Eqs. ([T2|) and 
(|g), this condition translates into A^e > 6(Tft™^^/10^^K)3/2(z/o/lGHz)3(i/b/1000z/o)°-^'^^ cm'^. 
Another possible cooling mechanism is adiabatic expansion of a plasma cloud. The charac- 
teristic time scale for this process is tad = 3 ■ 10^°(i?/lpc)(t//10%ms~^)-is, where R is the 
size of the cloud and U is the expansion velocity (assuming spherical expansion). 
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4 Effect of induced Compton scattering on the electron 
distribution 



So far we have assumed that the distribution of the electrons in momentum space remains 
relativistic Maxwelhan during their interaction with the low-frequency radiation. This should 
be the case if the plasma is dense enough that electron-electron collisions can quickly smooth 
out any arising deviations from a Maxwellian distribution. In order to find out whether this 
thermalization really takes place, one should compare the induced Compton heating time, 
^heat, given by Eq. (^), with the time scale on which electrons can relax to a Maxwellian 
distribution (Spitzer, 1978), te-e = 2.5 ■ W'^{lnA/4:0)-\kT/mc^f/'^N-^ s, where InA is the 
Coulomb logarithm. Taking T = T^^ (Eq. [^), we find that relaxation is efficient when 

- c^y f^y f^v"'%- < 1. (28) 
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Let us consider a few examples. At T™^^ = 10^^ K and uq = 1 GHz, Maxwellization of the 
electron spectrum occurs when ^ 10~^ cm~^. For T™^^ = 5 ■ 10^^ K and z/q = 10 GHz, 
the corresponding range is ^ 10^ cm~'^. It should be noted here that the establishment of 
the high-velocity tail of the Maxwellian distribution occurs on a larger time scale than t^-e 
given above. In the nonrelativistic limit, the appropriate characteristic time is proportional 
to {v/ (v))^ when v ^ (v), where v is the velocity of an electron and {v) is the typical thermal 
velocity of electrons (e.g., Krall, Trivelpiece, 1973). 

We now consider an extreme situation when tg-c ^ ^heat, i-e., no Maxwellization of electrons 
takes place due to collisions. To find out what momentum distribution, f{p), results in this 
case, we need to consider the diffusion of electrons in momentum space caused by induced 
and spontaneous (inverse) Compton scattering. The corresponding Fokker-Planck equation 
for the electron momentum distribution is (Illarionov, Kompaneets, 1976) 



9/ _ l_d_ 
dt dp 



p 



dp 



(29) 



where the diffusion coefficient 



D{f3) 



2(3' 



74/38(1 + /?')' 



^{y')GD{P,P')dp' 



(3 - (3^)1^^-^+2(3- 2(3' r-] + 2(f + - /?') 

y V 7 / 



(30) 



and the breaking force due to spontaneous scattering (Landau & Lifshits 1975) 



--crTe/37^ 



(31) 



The electron velocity is related to the electron momentum via 



P 
•jmc 



P 

mc 



1 + 



{-} 

Kmc/ 



-1/2 



(32) 
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The remaining quantities appearing in Eqs. 



and (pTI) are the same as in Eq. 



The cumbersome expression ([30| ) can be simphfied if we assume that {3 <^ 1. Series expansion 
can then be carried out, which is completely analogous to those written down in §3 for the 
induced Compton heating rate: 



(TtC 

^ivT Jo 



21 S^{de J dufdu' 



+ 



(33) 



The equilibrium momentum distribution is given by 



feqip) = const exp 



(34) 



In the nonrelativistic limit, D{p) = const and /sp(p) ~ P', hence F^q ~ exp{—p'^/2mkT^^) 
with Tgq given by Eq. ([TBI) . Therefore, the equilibrium distribution is Maxwellian in this 
limit, which was first shown by Zel'dovich and Levich (1970). 

We can find the first-order relativistic correction to the nonrelativistic equilibrium distri- 
bution from Eq. (0), by making use of the approximate expression ( ^3]) for the diffusion 
coefficient and transforming Eq. (^) to Fgp = — 4/3crTe(p/mc)[l + 0.5(p/mc)^ + ...]. The 
result is 



feqip) = const exp 



p 



*=4 




2Aj£(^eJdufdi/ 
50 J,^ elu-^du , 



(35) 



where is given by Eq. (|T6| 



In the case of synchrotron radiation with self-absorption at low frequencies (Eqs. |T2[, |T^), 
the equilibrium distribution is 



feqip) = const exp 



p^ 



2mkT^„ 



1 + 0.69 



(-)'■ 

\mcj 



(36) 



In Fig. ^, we have plotted the equilibrium distributions for two values of T™^^: 5 ■ 10^^ K 
and 5-10^^ K. The corresponding stationary temperatures, as estimated using nonrelativistic 
Eq. (16), are /cTeq = 95 keV and 950 keV. The exact result was obtained by numerical eval- 
uation of Eq. (133) using Eqs. (|5D|) and (PT|). Also plotted are the nonrelativistic Maxwellian 



distribution at temperature T^^ and the mildly-relativistic approximation given by Eq. (0). 
One can see that the right wing of the distribution is substantially suppressed compared 
to the nonrelativistic Maxwellian distribution. This takes place because of the decreasing 
diffusion coefficient and increasing breaking force with increasing p. Surprisingly enough, 
the approximate formula (^) provides an excellent fit even when kT^^ ^ mc^. 

Since the distribution ( p4D becomes a Maxwellian one in the limit kT^^ -C mc^ and assumes a 
quasi-Maxwellian shape in the mildly-relativistic regime, it is natural to characterize this type 
of distribution by some effective temperature. Teg. We define Tgfj to be the temperature of the 
relativistic Maxwellian distribution for which the mean electron energy, (7)mc^=((p^ + l)^/^), 
is equal to that for a given quasi-Maxwellian distribution, i.e.. 



I P^ feqip) dp 



(37) 
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where the dependence of (7) on temperature is given by Eq. (|T^). 

For the two example distributions presented in Fig. ^, kT^s takes values of 45 and 140 keV 
— these should be compared with kT^^ = 95 and 950 keV, respectively. The relativistic 
Maxwellian distributions that correspond to these Tgfj values are shown in Fig. ^. 

In Fig. ^, we have plotted T^s, calculated using both the exact formula (0) and its mildly- 
relativistic approximation (|35D , as a function of T™^^. For comparison, the dependence 
Teq(T™'^^) for Maxwellian electrons is reproduced from Fig. ^. One can see that the two exact 
dependencies are nearly coincident (the difference is less than 10%) in an extremely broad 
range of parameter values: fcTefj, fcTgq ^ mc^. Even more surprising is the nearly perfect 
agreement between the exact solution for Teq(T™'^^) and the midly-relativistic approxima- 
tion (Eq. m) for Tefr(T^°'^^). We should emph asize here that the differences between the 



different dependences shown in Fig. ^ are real (as is confirmed by the fact that they diverge 
significantly when Tes ^ mc^), although very small. An important conclusion follows: the 
energy that can be accumulated by an ensemble of electrons as a result of induced Compton 
heating (with inverse Compton scattering serving to cool the electrons) is almost indepen- 
dent in the sub-relativistic regime on whether the electrons are maintained Maxwellian while 
being heated or not. 

We have also checked that the contribution of the heated electrons to the gas pressure, which 
is proportional to (pP) ~ /(p^ + 1)^/^[1 — (1 + p'^)~^]p'^ f {p)dp (compare with Eq. for the 
mean energy), proves to be nearly the same (to within 2% for arbitrary values of Tes) for 
a quasi-Maxwellian plasma with the distribution function ( ^4|) and for the thermal plasma 
with T = Teff (note that in the limit kTes ^ mc^ both pressures must be equal, because 
j3 ^ 1). Therefore, the effective temperature Teg perfectly characterizes the thermodynamic 
properties of a mildly relativistic quasi-Maxwellian plasma that has been heated by means 
of the induced Compton process. 



The above discussion also suggests that one can make use of the simple expression (|35D to 
estimate the mean stationary electron energy with high accuracy in a very broad parameter 
range, at least up to Teg ~ mc^, which for our model spectrum corresponds to as high as 
jimax ^ ]^gi5 independent of whether thermalization takes place or not. 



Suppose that plasmas whose electrons are in the distribution given by Eq. (|35| ) do exist. 
Would they be different from thermal plasmas observationally? In particular, will the energy 
spectrum of hard X-ray bremsstrahlung emission from such plasmas be peculiar? We have 
computed such a spectrum for the momentum distribution shown in Fig. |a (/cTeg = 45 keV). 
The result is presented in Fig. ^ The computation consisted of weighting the Bethe-Heitler 
formula (Jauch, Rohrlich, 1976) for the cross-section of electron-ion bremsstrahlung with 
the given momentum distribution (Eq. |^). One can see that the deviation from the 
spectrum produced by an ensemble of relativistic Maxwellian electrons with temperature 
Tefi is negligible. Only when kT^^ ^ 100 keV, does the bremsstrahlung spectrum formed 
in the quasi-Maxwellian plasma become noticeably different from that corresponding to a 
thermal plasma. 

This research was partially supported by the Russian Foundation for Basic Research (projects 
00-02-16681 and 00-15-96649). 
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Appendix 



In a recent paper (Sazonov, Sunyaev, 2000) we have given an analytic expression for the 
kernel of the integral kinetic equation describing the redistribution of photons in frequency 
as a result of induced Compton scattering in a mildly relativistic thermal plasma. This 
kernel allows one, in principle, to derive the terms associated with the induced Compton 
process in Eq. (HD — the generalized Kompaneets equation. Unfortunately, the published 
expression contains a minor error (in the leading coefficient), and as a consequence, the 
formula does not correspond to the stated accuracy (it is valid when kT ^ O.Olmc^, instead 
of kT <^ O.lmc^). We use the opportunity to give here the correct expression: 



kT 



327r V vr \mc^ 
-/i(l - ^2) 



-3/2 



hv'{v' - v) 



mc^gu 



1 + /i^ + - - /i 



63 



kT 



3(1 + ij,^)mc^ v' -V 



8kT 







)*1 


exp 



2g^kT 



2vu' jJL + u 



/2\l/2 



(38) 



Here, v and are the frequency and direction of propagation of a photon before scatter, v' 
and Vt' are the corresponding values after scatter, and ^ = flQ,' is the scattering angle. 

An analogous correction is due for the kernel averaged over the scattering angle (see Sazonov, 
Sunyaev, 2000): 



P^°^(z/; u') 

Po 

Pt 
F 



u) -Piu 



V 



^ kT\ hu'(u' 



V 



m(?v{i> + v') 
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1091 507 
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68 



) G, 



exp 



1120 560 

(-6^), G 



35 



35 



26 



136 



35 



kT 



/ exp{-e)dt = 0.57r^^'^Erfc{\5\), 

J \6\ 



kT 



-1/2 , 



V 
V 



(39) 



The formulae (|38| ) and (p9D are applicable in the range hv <^kT ^ O.lmc^ nhv = kTb > kT 
(the latter condition means that the inverse Compton effect is small with respect to the 
induced one). 



References 

Bicknell G.F., Dopita M.A., O'Dea C.P.O.// Astrophys. J., 1997, 485, 112 
Blandford R.D.// Astron. Astrophys., 1973, v. 26, p. 161. 



14 



Blandford R.D., Scharlemann E.T./ / Astrophys. Space Science, 1975, v. 36, p. 303. 
Challinor A., Lasenby A.// Astrophys. J., 1998, v. 499, p. 1. 

Celotti A., Kuncic Z., Rees M.J., Wardle J.F.C.// Mon. Not. R. Astron. Soc, 1998, v. 293, 
p. 288. 

Gould R.J.// Astron. Astrophys, 1979, v. 76, p. 306. 

niarionov A.F., Kompaneets D.A.// Sov. Phys. JETP, 1976, v. 44, p. 930. 
Itoh N., Kohyama Y., Nozawa S.// Astrophys. J., 1998, v. 502, p. 7. 

Jauch J.M., Rohrhch F.// The Theory of Photons and Electrons 2nd ed.). New York: 
Springer- Verlag, 1976. 

Kellerman K.I., Vermeulen R.C, Zensus J. A., Cohen M.H.// Astrophys. J., 1998, v. 115, 
p. 1295. 

Kompaneets A.S.// Sov. Phys. JETP, 1957, v. 4, p. 730. 

Krall N.A., Trivelpiece A.W./ / Principles of plasma physics. New York: McGraw-Hill, 1973. 

Kuncic Z., Bicknell G.V., Dopita M.A.// Astrophys. J., 1998, v. 495, L35. 

Landau L.D., Lifshitz E.M.// The Classical Theory of Fields (4th ed.), Oxford: Pergamon 
press, 1975. 

Levich E.V., Sunyaev R.A.// Sov. Astron., 1971, v. 15, p. 363. 
Peyraud J.// J. de Phys., 1968, v. 29, p. 88. 

Pozdnyakov L.A., Sobol I.M., Sunyaev R.A.// Astrophys. Space Phys. Rev., 1983, v. 2, p. 
189. 

Sazonov S.Y., Sunyaev R.A.// Astrophys. J., 2000, v. 543, p. 28. 

Shestakov A.I., Kershaw D.S., Prasad M.K.// J. Quant. Spectr. Rad. Transf., 1988, v. 40, 
p. 577. 

Sincell M.W., Krolik J.H.// Astrophys. J., 1994, v. 430, p. 550. 

Spitzer L./ / Physical Processes in the Interstellar Medium, Chichester: Wiley, 1978. 

Sunyaev R.A.// Sov. Astron., 1971, v. 15, p. 190. 

Vermeulen R./ / in Proc. "Astrophysical Jets" , Cambridge: Cambridge Univ. Press (eds. 
Burgarella D. et al.), p. 241, 1992. 

Vinogradov A.V., Pustovalov V.V.// Sov. Phys. JETP, 1972, v. 35, p. 517. 

Woodward P.// Phys. Rev. D, 1970, v. 1, p. 2731. 

Zel'dovich Ya.B.// Sov. Phys. Usp., 1975, v. 18, p. 79. 

Zel'dovich Ya.B., Levich E.V.// Sov. Phys. JETP Lett., 1970, v. 11, p. 35. 



15 



0.1 



0.01 



0.001 



1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 




J I I I M 1 1 1 I I I I I I I I I I I I I 1 1 1 1 I I I 1 1 1 



0.1 



10 100 
Frequency 



1000 10000 



Figure 1: Radiation spectrum produced by a self-absorbed synchrotron source, described by 



Eqs. ( |12[) and (|13|) . The frequency is measured in units of the characteristic frequency uq. 
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Figure 2: (a) Deviation of the induced Compton heating rate of Maxwelhan electrons located 
in an isotropic field of self-absorbed synchrotron radiation from the nonrelativistic estimate, 
W(l (Eq. [[l^]), as a function of temperature. The solid line — the exact result, obtained by 
weighting Eq. (^) with a relativistic Maxwellian distribution. The dashed and dash-dotted 
lines represent the results of the calculation by the mildly-relativistic formula (|1^) in which 
retained are, respectively, only the correction term 0{kT /mc^) and all quoted terms up to 
0{{kT/mc^f). (b) Deviation of the inverse Compton cooling rate from the nonrelativistic 
estimate, Wq{T) (Eq. ||l9l). The solid line — the exact result (Eq. fl^l)- The dashed and 
dash-dotted lines represent the results of the calculation by the mildly-relativistic formula 



op, similarly as in (a). 
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Figure 3: Stationary electron temperature vs. peak radiation brightness temperature. This 
dependence resuhs from the balance W~^(Teq) = W~(Teq), with W^(T) and W~{T) as 
plotted in Fig. ^. The types of the lines have the same meaning as in Fig. |^. Also shown 
(the dotted line) is the nonrelativistic result (Eq. [|T6[ or Eq. 
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Figure 4: Evolution of the temperature of Maxwellian electrons during their Compton in- 
teraction with self-absorbed synchrotron radiation. The curves are labelled with the cor- 
responding values (in K) of the peak brightness temperature T™^^. At moment t = the 
plasma is cold. 
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Figure 5: (a) Equilibrium distribution of electrons located in an isotropic field of self- 
absorbed synchrotron radiation with T{^^^ = 5 ■ 10^^ K. The solid line — the exact result, 
obtained using Eqs. (P^), ( PO]) and (^TJ). The dotted line — the nonrelativistic Maxwellian 
distribution with fcT^q = 95 keV, a value found from Eq. ([T6|). The dashed line (almost 
coincident with the solid line) — the mildly-relativistic approximation (Eq. ||36|). The dash- 
dotted line — the relativistic Maxwellian distribution with kT = kT^s = 45 keV, a value 
found from Eq. (|37D . The mean electron energy for this distribution is equal to that for 
the equilibrium distribution, (b) The same as (a), but T™^^ = 5 ■ 10^^ K, in which case 
kT^ = 950 keV and kT^f, = 140 keV. 
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Figure 6: Effective temperature of tlie stationary electron distribution, defined by Eq. (^), 
vs. the peak radiation brightness temperature. The sohd hne — the exact resuh, obtained 
using Eqs. (|3^) , (pOD and (|3TD . The dashed hne — the result for the distribution given by 
the approximate formula (^). The dependence of the stationary temperature of Maxwellian 
electrons on T™^^ is reproduced from Fig. ^ for comparison (the dash-dotted line). 
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Figure 7: Spectrum of hard X-ray bremsstrahlung emission from a plasma whose electrons are 
in the distribution plotted in Fig. ^ (the solid line), compared with the spectrum (the dashed 
line) that corresponds to the relativistic Maxwellian distribution with kT = kT^s = 45 keV. 
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